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Abstract 

In the context of quasi-Hermitian theories, which are non-Hermitian 
in the conventional sense, but can be made Hermitian by the intro- 
duction of a dynamically-determined metric r], we address the prob- 
lem of how the functional integral and the Feynman diagrams de- 
duced therefrom "know" about the metric. Our investigation is trig- 
gered by a result of Bender, Chen and Milton, who calculated per- 
turbatively the one-point function Gi for the quantum Hamiltonian 
H = ^{p^ + x"^) + igx^. It turns out that this calculation indeed cor- 
responds to an expectation value in the ground state evaluated with 
the r] metric. The resolution of the problem turns out be that, al- 
though there is no explicit mention of the metric in the path integral 
or Feynman diagrams, their derivation is based fundamentally on the 
Heisenberg equations of motion, which only take their standard form 
when matrix elements are evaluated with the inclusion of r]. 

PACS numbers: 03.65. Ca, ll.30.Er, 02.30.Cj 

1 Introduction 



The original paper of Bender and Boettcher[T] gave convincing numerical and 
analytic evidence that a class of non-Hermitian Hamiltonians can nonetheless 
possess a completely real spectrum, in their case because of an unbroken PT 
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symmetry. However, at first sight such theories, treated as fundamental quan- 
tum mechanical theories, are still not acceptable because the natural Hilbert- 
space metric, whereby the overlap between two wave functions x ^i-iid ip is 
not the standard J dxx*{x)ilj{x), but rather J dxx*{—x)ilj{x) = f dx{x)pT'ip, 
is non-positive-definite and so does not lend itself to the usual probabilistic 
interpretation. Such a probabilistic interpretation was subsequently resur- 
rected by the realization [21 E] that a positive-definite metric could in fact be 
constructed, although it is not determined once and for all, like the standard 
metric, but depends on the Hamiltonian itself, and therefore differs from one 
model to another. 

The authors of Ref. [2] introduced the so-called C operator, with coordinate- 
space representation C{x,y) = Ylin't'n{.x)4'n{.y) in terms of the PT-invariant 
eigenstates (pn^x) of H, and the overlap / dx{x)cPT{.x)ip{x) , where {x)cpt{.x) = 
J dy C{x,y)x*{y)- It was subsequently [1] found useful to write C in expo- 
nential form, C = e^P. The Q operator provides the link with the notation 
of MostafazadelJll, namely i] = e~'^ . In this notation, the matrix element of 
any operator O is 

{{x\Om = {X\V0\^P). (1) 
C commutes with H, which translates into the statement of quasi- Hermiticity: 

ift = e~^He^ = 7]H7]-\ (2) 

As realized by Mostafazadeh, this latter equation can be used to define an 
equivalent Hermitian Hamiltonian h, with the same spectrum as H, by the 
similarity transformation 

h = e-^^He^^ = pHp'^. (3) 

The terminology "quasi-Hermitian" is now commonly used to describe non- 
Hermitian Hamiltonians H thus related to a Hermitian Hamiltonian h. The 
metric r] is central to the formulation of such theories in the Schrodinger 
formalism. However, its role in the functional integral formalism is much less 
prominent, even subliminal[Sl E]- Nonetheless the functional integral, and 
the Feynman rules derived from it, do somehow "know" about the metric. 

The present investigation is prompted by the calculation by Bender et 
al. [7] of the one-point function for the imaginary cubic Hamiltonian, 

H = ^{p^ + x^) + zgx\ (4) 

^Strictly speaking we should write 77+ for 77, to emphasize the fact that it is positive- 
definite 
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corresponding to the Lagrangian 



(5) 



using the standard Feynman rules derived therefrom. The result, to order 
is 



It turns out that this is the result of calculating the Green function, using 
the formula of Eq. ([T]), as the expectation value of = a; in the ground state 
\x) = IV') = l^o)- In other words, the graphical calculation of Gi corresponds 
specifically to the matrix element calculated using the t] metric, even though 
that metric appears nowhere in the former calculation. 

It is the purpose of the present paper to elucidate this problem, and 
answer the question posed in the title. The plan of the paper is as follows: in 
the next section we look at a simpler problem, the Swanson model [H], where 
we do not even have to deal with perturbation theory. In Section 3 we discuss 
in more general terms the relations between the Schwinger-Dyson equations 
and perturbation theory for the imaginary cubic potential. In Section 4 we 
see how the same results arise from a functional integral. In Section 5 we 
summarize what we have learned from this analysis. The Appendix deals 
with the relation between the Schrodinger and Heisenberg pictures for quasi- 
Hermitian theories, which turns out to be crucial. 

2 The Swanson Model 

The classical Hamiltonian for this model is 



where a and b are real and positive, while c is pure imaginary. The construc- 
tion of Q is not unique. 



We first suppose that Q = Q{x). The associated Lagrangian is obtained by 
setting X = OH /dp = 2{bp + cx), to give 




(6) 



H = ax^ + hp^ + 2cxp, 



(7) 



2.1 Q = Q{x) 



where a = a — c^/b, after having discarded a total time derivative. This latter 
is in fact a multiple of the derivative of the Q operator, when this is taken 
to be a function of x only [5]. Thus the Lagrangian L associated with H is 
the standard Hermitian Lagrangian of a scaled harmonic oscillator. Hence 
both Feynman diagrams and Green functions may be evaluated within the 
standard framework of quantum theory, for which the Hermitian Hamiltonian 
h corresponding to L is (with P = dL/dx = x/ {2b) = p + {c/b)x) 

h{x,P) =bP'^ + ax^. (9) 

Therefore the Green functions . . . tn) are just standard vacuum expec- 

tation values of time-ordered products in the vacuum state. For simplicity 
let us just consider G2(t,t), which is 

G2{t,t) = {Qhlx'lQh) (10) 

The ground state of the original non-Hermitian H is related to that of h by 
the similarity transformation of Eq. (jS]), namely 

\Qh) = e~'^'^\nH) . (11) 

Thus 

G2(t,t) = (l^H|e-^«xV^«|fiH) (12) 

In this case, because Q = Q{x) it commutes with x, so that 

G2{t,t) = (finle-^x^ll^n) = mulx'inH)) (13) 

Thus the Green function associated with the Lagrangian L is unequivocally 
the vacuum expectation value of x evaluated using the ry metric. 

2.2 Q = Q{p) 

Although Q = Q{x) was the simplest possibility, there is a one-parameter 
family of Q operators P, [lO], where Q can depend on both x and p. Perhaps 
the extreme opposite of the previous case is when Q = Q{p). As we have 
seen, a straightforward calculation of the Lagrangian leads us to Q = Q{x), 
so how can we access the other cases? The answer is that we must first make 
a change of variables. 

For the case Q = Q{p) we consider the Hamiltonian H[C,,p), with yet 
to be specified, and make the shift 

^ = X - {c/a)p . (14) 
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Then H becomes 

H{^,p) = ax^ + hp^ = h{x,p), (15) 
where b = b ~ c^/a, with corresponding Lagrangian 



^2 

X 2 



i(x) = ^-ax' (16) 

Now we have a new Hermitian harmonic oscillator, which, however, is a scaled 
version of the previous one, since ab = ab. Thus the ground states \Qh) are 
identical. 

Proceeding as in §2.11 to derive Green functions derived from the func- 
tional integral or Feynman diagrams for H{C,,p), we will obtain standard 
expectation values of products of ^(tr) with respect to the ground state \ flh), 
so in particular 

G2(t,t) = Wei^h) 

= (fiH|e-^«eV^«|fiH). (17) 

At this point it is useful to supplement Eqs. and with the more general 
statements that an observable A, with real eigenvalues, is quasi-Hermitian: 

= e-^Ae^ (18) 

and is related to its Hermitian counterpart a by 

A = e^/^ae-Q/^, (19) 

or equivalently 

At = e-«/2aeQ/2. (20) 

In the present case Q{p) is actually Q = i{c/a)p'^, which gives X = x+{c/a)p 
for the observable associated with x. Thus we see that in Eq. 0141) the variable 
^ must be identified with X"^ . Accordingly 

G2{t,t) = {QH\e~'^{e-^^{X^fe"-^'^)\nH) 

= {nH\e"'^x^\nH) = {{nHlx^inn)) . (21) 

What we have effectively done is to make use of the second of the identities 
H{x,p) = e-^/^h{x,p)e^/^ = h{X,P), 

(22) 

H{X\P^) = e~^/^h{X\P^)e^/^ = h{x,p), 
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where, in this particular case, P = p. A completely parallel treatment of 
case 12.11 would have started with H{x,tt), where tt = P^. However, we did 
not do this because the momenta drop out of the calculation immediately L 
is formed. With this understanding of the rj metric in mind, we return to the 
problem as posed, the Feynman diagram representation of Green functions. 

3 The Schwinger-Dyson Equations 

We consider a general non-Hermitian Hamiltonian of the form 

H=^{p'' + x^) + Uix). (23) 

Let us assume for the moment that the Green functions are calculated using 
the 7] metric, and verify later that that is indeed the case. Then the vacuum- 
generating functional is 

Z[j] = {n\r]T{exp[z J dtj{t)xm\n), (24) 

which by successive functional differentiation with respect to the external 
current j(t) gives 

Gniht2...tn) = {n\7]T{xiti)xit2)...xitn))\n) 

(25) 

1 



m 



-t5/5jit,))i-^S/63it2))...i-^6/5jit^))Z[j] 



j=0 



For a given U the Schwinger-Dyson equations are derived from the operator- 
valued equation of motion (see App. A) 

{d'^ + l)x{t) + U'{x{t)) = 0. (26) 

Since rj is time independent this gives rise to a functional differential 
equation for Z[j], namely 

[{d^ + l)H6/6j{t)) + U'{-z6/6m)] Z[j] = mZ[j]. (27) 

We stress that this last equation contains no reference to r] and is just as 
we would expect for a standard Hermitian theory. However, had we not in- 
cluded the correct metric r] we would not have obtained the usual Heisenberg 
equations of motion. 
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For the case of interest, U = igx^, we can obtain the Schwinger- Dyson 
equations by repeated functional differentiation of flTH) with respect to at 
j = 0. In terms of the propagator A(t — t'), the solution to {d^ + l)A(t — t') = 
—6(t — t') with appropriate boundary conditions, the first few are: 

= 3ig j dt'A{t-t')G2it't') = 0, 
G2{tti) = iA{t-ti) + ?>ig j dt'A{t-t')G-i{t't'ti), (28) 
G^{ttit2)) = i[A{t - ti)Gi{t2) + (ti ^ ts)] + "iig j dt'A{t - t')Gi{t't'tit2), 

Exphcitly, A{t) = -|ze-*l*l. Note that Gi{t) = {{x{t))) = x is in fact a 
constant by time-translation invariance. 

Here the Gn are not connected. The relations to the connected Green func- 
tions Wn, with generating functional 

W[j] = -z\nZ[j], (29) 

are, with Gi = Wi = x, and 14^2 («') = W2it- t'), 

x = 3ig{j dt'A{t - t'))[x^ + 1^2(0)] = ~3ig[x'^ + ^2(0)], 

W2{t - ti) = iA{t - ti) + 3ig j dt'A{t - t')[2xW2{t' - h) + ^3(^1)], (30) 

W3{ttit2)) = 3ig / dt'A{t - t')[2xW3{t'tit2) + 2W2{t'ti)W2{t't2) + W4{t't'tit2)], 



In the final equation we expanded about the time t. All these, and further 
equations, have simple diagrammatic representations. We stress that these 
are the conventional combinatoric equations for connected Green functions 
in an x^ theory with imaginary coupling constant. 

The Schwinger-Dyson equations (l30l) do not by themselves give a unique 
solution. Thus, a priori Wi is completely undetermined, and W2 is only par- 
tially. However, they have a unique perturbative expansion in g [TT]. The 
lowest terms in standard perturbation theory follow from the first and second 
of Eqs. ( l30i) . namely 

W2{t) = zA(t)+0(/) 

X = -3zgW2{0) + O{g') = -^zg + O{g^), (31) 
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the latter corresponding to the tadpole diagram of Ref . [7] . 

As mentioned in the Introduction, and as already remarked in Ref. [H], 
this value unambiguously corresponds in the Schrodinger picture to the ex- 
pectation value 

= Jdx - gQ,)xMx), (32) 

where 

Qi = -^p^-2xpx (33) 

12 1 

and ipo{x) = e~2^ /ni. Its non-zero value is entirely due to the presence of 
the e~'^ factor. We have also checked the 0{g^) result reported in Eq. dU]). In 
that case we need both Q3, occurring in the expansion Q = gQi+ g^Qs + . . . , 
and the third-order wave function. 

We can now identify the reason why the perturbation expansion of the 
Schwinger-Dyson equations and the resulting Feynman diagrams correspond 
to matrix elements calculated using the 77 metric. As detailed in Appendix A, 
it is only when matrix elements are so calculated that the Heisenberg-picture 
fields obey the standard equations of motion, which are the starting point 
for the Schwinger-Dyson equations and, subsequently, the functional-integral 
representation of Green functions and the Feynman diagram expansion. 

Note that as an expectation value in the ground state \(po) of the equiv- 
alent Hermitian Hamiltonian 

X = {^po\e~'^x\ipQ) 

= {(po\e^^e~^xe^'^\(fo) 

= (<^o|^Vo) (34) 

This latter calculation was also performed diagrammatically in Ref. [7j. Their 
variable x is in fact 

In the present case there is another, non-perturbative, solution to Eqs. fl3UI) 
starting with x = i/{3g), which is actually the lower minimum of the classical 
potential. However, after making the shift 

i 

x = — + y (35) 
3^ 

the Hamiltonian becomes 

H = ^{p^- y^) + igy"" + const., (36) 

which corresponds to a theory with the wrong sign of the quadratic term and 
is hence unacceptable. 
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4 The Functional Integral 



The formal solution to Eq. ( 1271) in terms of a path integral is 



= / [Dx] exp 



dt 



1.2 1 2 

■X X 



U (x) + jx 



(37) 



2 2 

In general one might have to worry about possible contours in the complex 
X plane on which the functional integral is defined. However, this is not 
an issue if we are restricting our attention to perturbation theory, since by 
Eq. ( l37j) we mean the perturbative expansion 



m = J [Dx] J2 



-t/(x))" 



exp 



dt 



1.2 1 2 

-X X 

2 2 



JX 



(3^ 



with each term in the series a derivative of the well-defined Gaussian integral 



Zo[j] 



[Dx] exp <i dt 



-X 



-X 



exp 



1 

'-i 
2 



2 2 

dtdt'j{t')A{t-t')j{t) 



JX 



(39) 



Although fl58l) recreates the perturbation series for the Schwinger-Dyson 
equations trivially, there is still no guarantee that such a divergent series 
is asymptotic. For the case in point U (x) = igx^ we know that the metasta- 
bility of the corresponding real potential U{x) = gx^ enforces a lack of 
Borel summability [I2] but we expect the complex potential to give Borel 
resummable series. 

The above concerns expectation values of products of x(tj). However, 
suppose we wished to do the same for other variables, in particular the cor- 
responding observable X, which typically is a function of both x and p. In 
that case we would have to start with the equivalent Hermitian Hamiltonian 
h{x,p), and exploit identities such as 



{{Qh[X^[Qh)) = Wx^[^h) 



(40) 



Because h{x,p) typically contains higher powers of p, one needs to use the 
full phase-space functional integral representation of Z[j], namely 



Z[j] = / [Da;][Z)p] exp < —i dt [px + h(x,p) — j x] 



(41) 



subject to the usual ordering problems. To rewrite this in terms of H one then 
needs to change variables to ^ = X\ vr = P"!" using the identity H{C,, vr) = 
h{x,p). This is what was done for the Swanson model in Ref. [5]. In Section 
2 we were instead concerned with expectation values of products of x and so 
started with {Qh[^'^[^h) ■ 
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5 Discussion 



Our main purpose in this paper was to understand how the perturbative 
calculation by Bender et al. [7j picked out a particular version of the one- 
point Green function for the non- Hermit ian theory with U = igx^, namely 
Gi = {ipolrjxlipQ) , rather than (ipolxlipo) , in spite of the fact that the metric 
rj = does not appear anywhere in the Feynman diagram expansion. 

After looking at the simpler soluble problem of the Swanson Hamiltonian, 
where we verified that the same feature occurred there for the two-point 
function G2{x,x), we went back to the derivation of the Schwinger- Dyson 
equations, which in turn lead to the functional-integral representation of 
the vacuum generating functional Z[j] and W[j], and to their perturbative 
expansions. Even then, the solution is not necessarily unique because of the 
problem of giving meaning to a divergent series (see [H] for a more detailed 
discussion) . 

Starting from the Heisenberg equation of motion for the field x{t), Eq. fpU|) . 
it is in fact the case that one would obtain the functional differential equa- 
tion for Z[j], Eq. ( |2?|) . whether or not the metric i] was included in the 
definition of Z[j] in Eq. ( 124|) . The crucial point is rather that, as shown 
in the Appendix, these equations of motion will only take their standard 
form iAfj = [Ah,H] in a quasi-Hermitian theory, if the relation between the 
Schrodinger and Heisenberg pictures is built on matrix elements defined with 
the inclusion of the metric. 

Whether or not these Green functions, being time-ordered expectation 
values of products of the non-observable field x{t), are useful is another 
question. It is worth remarking that in our previous paper [5] on functional 
integrals for non-Hermitian Hamiltonians we were instead concerned with 
the functional integral representation of Green functions of observables, such 
as {{Qh\X'^\^h)) in Case 2.2 for the Swanson Hamiltonian. 

Finally, there is nothing peculiarly quantum mechanical, as distinct from 
quantum field theoretical, about our analysis, which extends naturally to 
pseudo-Hermitian quantum field theory. 

Appendix 

Schrodinger and Heisenberg Pictures 

Let us denote a Schrodinger-picture state by \ip) and its corresponding 
Heisenberg state by so that 




iHt 




(42) 
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identifying the two pictures at t = 0. Correspondingly 

{^\ = (^|e^^^* (43) 

We identify the Heisenberg operator Anit) corresponding to a given Schrodinger 
operator A by considering 

{^p\vA\ip) = {^\e'^^'r]Ae-'^'\^) 

= (y.|r/e^^*Ae-^^V), (44) 

by virtue of Eq. (j2]). Thus we identify 

Ah = e'"'Ae-'^\ (45) 

the usual transformation. However, it is no longer unitary because H is not 
Hermitian. Since we have no explicit time dependence in the theory, rj itself 
is time independent because Q^H'ft^^-iHt _ From Eq. fl45p we obtain the 
Heisenberg equation of motion for Ah'- 

tj^AH = [AH,H]. (46) 
This equation does not, however, apply to rj. 
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